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Abstract 

We update our search for the lepton flavor violating t~ /^~7 and r~ e~7 decays based on 
535 fb~^ of data accumulated at the Belle experiment. No signal is found and we set preliminary 
90% confidence level upper limits: B{t- ^"7) < 4.5 x 10"^ and B{t- ^ e-7) < 1.2 x lO"'^. 

PACS numbers: 13.35.Dx, ll.30.Fs, 14.60.Fg 
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INTRODUCTION 



To search for new physics beyond the Standard Model, we have been looking for the 
lepton flavor violating (LFV) r ^ fi'-f and r ^ 67 decays in the Belle experiment [1] at 
the KEKB asymmetric-energy e~^e~ collider [2]. Previously, we obtained the upper limits of 
B{t ^ fi-f) < 3.1 X 10"^ [3] and B{r 67) < 3.9 x 10"^ [4] at the 90% confidence level (CL), 
using about 86 fb~^of data recorded at the T(45') resonance. Later, the BaBar collaboration 
obtained the upper limits B{t < 6.8 x 10~^ [5] and B{t 67) < 1.1 x 10"''' [6] with 

232.2 fb~^of data. Here we report our updated analysis with a data sample of 535 fb~^. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [1]. 

The basic analysis procedure is similar to our previous one [3, 4]. The selection crite- 
ria for r i^ll&l are determined and optimized by examing Monte Carlo (MC) simu- 
lated singal and background (BG) events, including the generic r"'"r~, qq (g = u,d, s,c,b), 
Bhabha, fj^'^fi~, and two-photon events. The BG r+r" events are generated by the 
KKMC/TAUOLA [7] and the response of the Belle detector is simulated by the GEANT3 [8] 
based program. 

Photon candidates are defined as in Ref. [9]. Muon candidates are identified by using 
a muon likelihood ratio, £^ [10], which is based on the difference between the range of 
the track calculated from the particle momentum and that measured by the KLM, which 
includes the value of formed from the KLM hit locations with respect to the extrapolated 
track. Identification of electrons is performed using an electron likelihood ratio, He, based 
on the dE/dx information from the CDC, the ratio of the energy deposited in the ECL to 
the momentum measured by both the CDC and the SVD, the shower shape in the ECL, the 
hit information from the ACC and the matching between the positions of the charged track 
and the ECL cluster [11]. 

Event Selection 

We select events that include exactly two oppositely-charged tracks and at least one 
photon, consistent with r+r~ decays: one t"^ (signal side) decays to ^'^'^ and the other (tag 
side) decays to a charged particle that is not a muon (denoted hereafter as i/i), neutrino(s) 
and any number of photons. 

Each track should have a momentum p'-^^ < 4.5 GeV/c in the center-of-mass (CM) frame 
and transverse component to the beam axis pt > 0.1 GeV/c within the detector fiducial 
region —0.866 < cos^ < 0.956 to avoid contamination from Bhabha and fi~^fi~ events. 
(Hereafter all the variables defined in the CM frame have superscripts "CM.") Each photon 
is required to have an energy > 0.1 GeV within the fiducial region. The total energy in 
the CM frame should be E™1^ < 10.5 GeV to further suppress Bhabha and fi'^fi events. 
The magnitude of the thrust vector, constructed from all selected charged tracks and photons 
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above, is required to be in the range from 0.90 to 0.98 in order to suppress fi~^fi~ and qq 
events. (Fig. 1(a)). 

For muon identification, we require a likelihood ratio £^ > 0.95 and p > 1.0 GeV/c. 
On the tag side, a track with < 0.8 is defined as /i. The photon that forms a (/X7) 
candidate is required to have > 0.5 GeV and —0.602 < cos6'^ < 0.829 to remove spurious 
combinations of 7's. 

The opening angle between the /i and 7 of (fj.'j), 0.4 < cos^^^ < 0.8, is particularly- 
powerful in rejecting t~^t~ background, which contains vr'^'s from r decays. The sum of the 
energies of the two charged tracks and the photon of the (/i7), E^^, should be less than 9.0 
GeV to reject fi~^fi~ events. The opening angle between the two tracks should be greater 
than 90° in the CM frame while the opening angle between the fi and the boost direction 
of its mother r from the CM frame in the r rest frame is required to satisfy cos O^j. < 0.4, 
to remove combinations of /i and 7 from BG. 

The following constraints on the momentum and the polar angle of the missing particle 
are imposed: Pmiss > 0.4 GeV/c and —0.866 < cos^miss < 0.956. Here, Pmiss is calculated 
by subtracting the momentum of all charged tracks and photons from the beam momenta. 
To remove the t~^t~ background, a requirement on an opening between the tag-side track 
and the missing particle is applied, 0.4 < cos9^^_^ < 0.98. We calculate the missing mass 



squared on the tag side, ml = {E^^ — -Etag)^ — Pmiss' "where E^^ (-Etap-) is the sum of the 



energy of the signal side /x and 7 (the sum of the energy of all tag-side particles assuming 
rriTT for the tag-side track), and then require —1.0 (GeV^/c'^) < ml < 2.0 (GeV^/c^), as 
shown in Fig. 1(b). Finally, a condition is imposed on the relation between Pmiss and the 




0.9 
thrust 



FIG. 1: (a) Length of the thrust vector and (b) distributions for r — > fij. Dots are data, 
open boxes show the BG Monte Carlo(MC) distribution and shaded histograms are the signal MC. 
Arrows indicate the selected region. 

mass squared of a missing particle, m^j^g = -Emiss ~ Pmiss- Pmiss > — 5 x m^j^g — 1 and 
Pmiss > 1.5 X m^jgg — 1, where -Emiss is the sum of the beam energies minus the sum of all 
visible energy and is calculated assuming the muon (pion) mass for the charged track on the 
signal (tag) side, Pmiss is in GeV/c and mmiss is in GeV/c^. 
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Background contribution 

After the selection requirements described in the previous subsection, the dominant BG 
source is r+r~ events with the decay — > ^^v^Vj. and a photon from initial state radiation 
or beam BG. Other sources are the radiative process and t~^t~ events with — vr^z/^. 

Two variables are used to identify the signals: Minv, the invariant mass of (/i7), and 
AE = Efj_^ — -Ebeam' energy difference between the (/X7) energy and the beam energy in 
the CM frame, where the signal should have Minv ~ and AE ~ 0. The resolutions in 
Minv and AE are estimated by fitting asymmetric Gaussians to the signal MC distributions 
giving aJif/^""^ = 14.49±0.10/24.24±0.13 MeV/c^ and a^2P^°^ = 35.29±0.49/81.41±0.94 
MeV, where 

^high/low ^g^^g 

the standard deviation at the higher/lower side of the peak. 
To compare the data and MC simulation, we examine a ±5(j region with 1.65 GeV/c^ 
< Minv < 1.85 GeV/c2 and -0.41 GeV < AE < 0.17 GeV, as seen in Fig. 2(a). A blind 
analysis method is taken: a slanting ±3cr ellipse region, indicated by the dashed ellipse in 
Fig. 2(a), whose detection efficiency is e^a = 6.05%, is covered till the final stage of the 
analysis. 

After the selections, we find 71 events remaining in data and 73.4±6.7 events remaining 
in MC in the ±5a region outside of the blinded ellipse. The latter is dominated by r+r~ 
events with initial state radiation, r+r~7, and is comprised of 58.8 ±4.3 (70.3 ±4.7) r~'"r~7 
events, 13.1 ± 4.9 (15.0 ± 5.3) /i+yU~7 events, with incorrect identification, and 1.6 ± 1.6 
(3.2 ±2.2) two-photon events, where the numbers in the parentheses are the BGs that remain 
in the entire sample. 

This background composition was understood in the previous analysis; the t^t^^ process 
yields contributions in the AE < region, while fi'^fi~'y events mostly have AE > 0. The 
latter rate is ~20% of the former, and one additional small contribution of ~5% is known to 
exist. This BG distribution is well represented by a combination of Landau and Gaussian 
functions, as found in Ref. [4]. The final event distribution in this data sample is found to 
follow well the BG function; the /i"'"/i~7 contribution is (20±13)% of the r"'~r~7 contribution, 
while the rest ((5.2 ±4.1)%) is from the two-photon process e+e~ e~^e~fi^fi~. 

Signal extraction 

After unblinding, we find 23 and 94 data events in the blinded and ±5a regions, respec- 
tively, while 15.0 ± 3.1 and 88.4 ± 7.4 events are expected from the MC. Figure 2(a) shows 
the event distributions in the Miy^^-AE plane. 

In order to extract the number of signal events, we employ an unbinned extended maxi- 
mum likelihood (UEML) fit with the following likelihood function: 

-(s+b) N 

^ = ^FT n i-'S, + bB,) . (1) 
■ i=l 

Here, is the number of observed events, s and b are the numbers of signal and BG 
events to be extracted, respectively. Si and Bi are the signal and BG probability density 
functions (PDF), where i indicates the i-th event. Si is obtained from the signal MC, and 
Bi is the PDF for background mentioned above, whose distribution is concentrated around 
AE ~ —0.2 GeV, as indicated by the solid line in Fig. 2(a). To enhance the signal detection 
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FIG. 2: (a) Minv ^ distribution in the search for r fij in the itScr region. Dots are the data 
and shaded boxes indicate the signal MC. The dashed ehipse shows the 3a bhnded region and the 
dot-dashed ehipse is the 2a signal region. Dashed lines indicate the ±2(7 band of the shorter ellipse 
axis, projected onto the longer ellipse axis. The solid line indicates the dense BG region, (b) Data 
distribution within the ±2(7 band. Points with error bars are the data and the shaded histogram 
is the signal MC assuming a branching ratio of 5 x 10~^. The solid curve shows the best fit. 

sensitivity and to avoid this dense BG region, we use a ±2a ellipse as the signal region for 
the UEML fit. The result of the fit is s = -3.9, b = 13.9 with N = 10. 

Figure 2(b) shows the event distribution within the ±2(T band of the shorter ellipse axis, 
projected onto the longer ellipse axis, and the best fit curve. No events are found around 
the peak of the signal distribution. The negative s value is consistent with no signal. 

We examine the probability for obtaining this result and evaluate the 90% CL upper 
limit using a toy MC simulation. The toy MC generates signal and BG events according to 
their PDFs fixing the expected number of BG events [b) ai b = b while varying the number 
of signal events (5). For every assumed s, 10,000 samples are generated following Poisson 
statistics with means s and b for the signal and BG, respectively; the signal yield (s'^*") is 
evaluated by the UEML fit. To obtain the upper limit at 90% CL (§90) we take the s value 
that gives a 90% probability of to be larger than s. The probability to obtain s < —3.9 
is 25% in a case of a null true signal. In other words, due to BG fiuctuations a negative s 
value is possible with a large probability, although the physical signal rate is positive [12]. 

The toy MC provides an upper hmit on the signal at the 90% CL as Sqq = 2.0 events 
from the result of the UEML fit. We then obtain the upper limit on the branching fraction 
B9o{t ^7) at the 90% CL as 

S9o(r-/x7) = ^<4.1xlO-«, (2) 

where the number of r pairs produced is Nrr = 4.77 x 10^, and the detection efficiency for 
the ±2(j ellipse region is e = 5.07%. 

The systematic uncertainties for the BG PDF shape increase sgo to 2.2. The systematic 
uncertainties for e arise from the track reconstruction efficiency (2.0%), the photon recon- 
struction efficiency (2.0%), the selection criteria (2.2%), the luminosity (1.4%), the trigger 
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efficiency (0.9%), and the MC statistics (0.3%). All errors are added in quadrature to yield 
the total uncertainty of 4.0%. This uncertainty increases the upper limit of the branching 
ratio by 0.2% [13]. Since the angular distribution for r — > yU,7 depends on the LFV inter- 
action structure, we evaluate its effect on the result by assuming the maximum possible 
variation, V ± A interactions, rather than the uniform distribution so far assumed in the 
MC analysis. No appreciable effect is found for the upper limit. 

Finally, the following upper limit on the branching ratio is obtained: 

B{t /i7) < 4.5 X 10"^ at the 90% CL. (3) 



T — > 67 

For r — i> 67 we use a procedure similar to that for r fi'j. 
Event Selection 

We examine a t~^t~ sample, in which one r goes to an electron and a photon, and the 
other r decays to a charged particle, but not an electron (^), neutrino(s) and any number 
of photons. The selection criteria are basically the same as those for r /i7, so below we 
describe only the differences. 

An obvious difference is the replacement of the by an e on the signal side, and using a 
^ veto rather than a /i veto on the tag side. The electron on the signal side (57) is required 
to have Cc > 0.90 and a momentum p > 1.0 GeV/c, while the ^ on the tag side should 
have £e < 0.1. Minor differences in the kinematical selection include requirements on the 
missing mass squared on the tag side and the opening angle between the tag-side track 
and the missing particle on the tag side: —0.5 (GeV^/c^) < ml < 2.0 (GeV^/c^), and 0.4 
< cos6'^^g_^ < 0.99. The other requirements are the same as those for r ji-^. 

The Minv and resolutions are evaluated as at^^l^"^ = 14.76 ± 0.18/25.38 ± 0.38 

MeV/c2 and a^2P^°'" = 35.66 ± 0.62/89.98 ± 1.72 MeV, and then the ±5a region over 1.65 
GeV/c2 < Minv < 1.85 GeV/c^ and -0.45 GeV < AE < 0.18 GeV is assigned for the signal 
evaluation. A slanting ±3cr ellipse is also blinded. 

After the selection, we find 34 and 29.9 ± 2.8 events remaining in data and MC, respec- 
tively, in the ±5cr region, but excluding the blind one. The BG is from r+r~7 events; no 
other BG process is found to contribute. 

Signal extraction 

After opening the blind we find 13 and 55 data events in the blind and ±5cr region, 
respectively, while the MC predicts 8.1 ± 1.6 and 42.8 ± 3.7 events. Figure 3(a) shows the 
event distribution in the M\-a^-AE plane. 

The signal extraction process is the same as that for r ^ /i7, described in the former 
section. The BG is composed of (18 ± 18)% e"'"e~7 (radiative Bhabha), while the remainder 
is r"^r~7. No other background source is found in MC. The UEML fit over the ±2(T ellipse 
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FIG. 3: (a) M[^^ - AE distribution in the search for r — > 67 in the ib5cr region. Dots are the data 
and shaded boxes indicate the signal MC. The dashed elhpse shows the 3(T bhnded region and the 
dot-dashed elhpse is the 2a signal region. The dashed lines indicate the ±2(7 band of the shorter 
ellipse axis, projected onto the longer ellipse axis. The solid line indicates the dense BG region, 
(b) Data distribution within the ±2(7 band. Points with error bars are the data and the shaded 
histogram is the signal MG assuming a branching ratio of 5 x 10^^. The solid curve shows the best 
fit. 



region results in s = —0.14, h = 5.14 with = 5. The toy MC gives a probability of 48% 
for obtaining s < —0.14 in the case of a null signal. Figure 3(b) is the same as Fig. 2(b), 
but for the r ^ 67 case. The upper limit of §90 = 3.3 is obtained by the toy MC in the case 
of the UEML fit result. The upper limit on the branching fraction is calculated as 

B9o(r-e7) = ^<11.7xlO-«, (4) 

where the detection efficiency for the ±2a ellipse region is e = 2.99%. 

The systematic uncertainties are essentially the same as those for r —> fi'-f. minor differ- 
ences are the selection criteria (2.5%), and the trigger efficiency (2.0%). The total uncer- 
tainty amounts to 4.5%, and it increases the upper limit of the branching ratio by 0.2%. In 
addition, the systematic uncertainties for the BG PDF shape increase §90 to 3.4. Taking 
into account this systematic error, we obtain the 90% CL upper limit, 

B{t 67) < 12.0 X 10"^ (5) 



SUMMARY 

We updated our LFV searches for t fi'y and r — * 67 decays, based on 535 fb~^ of data, 
i.e. with about six times higher statistics than before: The resulting upper limits on the 
branching fractions are 

B{t ii-f) < 4.5 X 10~^ (6) 
B{t 67) < 12.0 X 10"^ (7) 
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at the 90% CL. 



For the t ^ fi'j search, we obtain e = 5.07% and A'^obs = 10 in the region where the 
upper hmit is evaluated, while we had e = 12.0% and A^obs = 54 in our previous analysis 
with the 86 fb~^ data sample. Similarly, for the r ^ 67 search, we obtain e = 2.99% and 
-^obs = 5 compared to e = 6.39% and A^'obs = 20 previously. Both modes have about 6 times 
better sensitivity than the previous analyses. We have also estimated the improvement in 
our sensitivity for the old data sample. The result would be only 2.5 times better if we 
applied the methods used in our previous analysis to the current data sample. 
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